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Abstract 

We present a program for calculating the isospin symmetry breaking of the B K*^ 
decay in the MSSM with minimal flavor violation. This program calculates the NLO 
super symmetric contributions to the isospin asymmetry, using the effective Hamiltonian 
approach and within the QCD factorization method. We show that isospin symmetry 
breaking proves to be a very restrictive observable, in particular in the mSUGRA param- 
eter space. The program also calculates the inclusive branching ratio associated to 6 — > 57 
transition, as a comparison reference. 

PACS numbers: ll.30.Pb, 12.15.Mm, 12.60.Jv, 13.20.He 



1 Introduction 

In spite of the numerous phenomenological achievements of the Standard Model, there are still 
theoretical motivations which impel particle physicists to believe in new physics beyond the 
standard model. Among these motivations, one can name in particular the hierarchy prob- 
lem or dark matter /energy problems. As an attempt to solve these problems, supersymmetry 
(SUSY) appears as a highly motivated and studied model during the past decades. 

We consider here the Minimal Supersymmetric extention of the Standard Model (MSSM), 
with minimal flavor violation. Since no supersymmetric partner has been discovered, super- 
symmetry has to be a broken symmetry. A large effort has been done to understand the nature 
of supersymmetry breaking, and the most successful mechanisms are believed to be the gravity 
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mediated (SUGRA) in which the SUSY breaking sector communicates with the MSSM sector 
via only gravitational interactions, anomaly mediated (AMSB) in which SUSY breaking sector 
is transmitted to the observable sector through the super- Weyl anomaly, and gauge mediated 
(GMSB) which is mediated by gauge interactions. The boundary conditions at high scales in 
the mentioned scenarios reduce the number of free parameters of the MSSM, allowing feasible 
phenomenological studies. 

During the past few years, rare B decays, and in particular those associated to 6 57 transi- 
tions have been extensively used in order to constrain SUSY parameter space. The observable 
commonly used for this purpose is the inclusive branching ratio corresponding to this transi- 
tion. In [1] we have introduced a new observable, the isospin asymmetry, and showed that this 
observable can produce complementary information to the branching ratio, and in some regions 
of the parameter spaces, even more restrictive constraints. 

We introduce here Super I so as a tool to calculate the isospin symmetry breaking for the three 
main SUSY breaking scenarios, i.e. mSUGRA, AMSB and GMSB. The program is modular 
so that the user can insert its own model (as long as it produces the SUSY mass spectra and 
couplings according to the SUSY Les Houches Accord conventions). The code allows also to cal- 
culate the inclusive branching ratio, so that comparison between the two observables is possible. 

The paper is organized as follows. The relevent equations for the calculation of the isospin 
asymmetry are summarized in section [2l In section [3] the Superlso package is described, and 
the installation and compilation instructions can be found in section HI In section [5] we explain 
how to use Superlso, with input and output samples. Finally some results are presented in 
section [HI 



2 Calculation of the isospin symmetry breaking 

The isospin asymmetry Aq- can be written to leading order [2]: 
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The functions Ki and can be written in function of the Wilson coefficients Q at scale fij^ 
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where N = 3 and Cp = 4/3 are color factors, and: 
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Here the number of flavors nf = 5, and fiQ = 0{mb) is an arbitrary normalization scale. 
The details of the calculation of Wilson coefficients is given in [1]. G^{xcb), H±{xcb) 

and X± are convolution integrals of hard-scattering kernels with the meson distribution am- 
plitudes, and their values are summarized in Table [TJ In this table also appears the param- 
eter X = ln(mB/A/i) (1 + ge^'^), which parametrizes the logarithmically divergent integral 
dx/ {1 — x). g < I and the phase if are arbitrary. Ah ~ 0.5 GeV is a typical hadronic scale. 



The coefficient reads 
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tpicase note that we adopt the conventions of [HIS], in which cf\Mw) = 1 and c'^\M\y) — 0, therefore 
Ci and C2 are inverted here compared to yy. 
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CKM parameters and B meson mass 


Vus 




Kb/Kb 


Re(lC.K6/K*.K6) 


TUB 


0.22 


0.041 ± 0.05 


0.085 ± 0.025 


0.011 ±0.005 


5.28 GeV 



B meson parameters 


fB 


\b 




< >_L 




200 ± 20 MeV 


350 ± 150 MeV 


0.19 ±0.02 


3.7 ±0.04 


(4.8 ±0.5) 
± (1.5±0.2)i 


K* meson parameters 


Fk* 


Ik* 


mK* 


Ik* 




0.38 ± 0.06 


175 ± 9 MeV 


892 MeV 


226 ± 28 MeV 


0.30 ±0.05 



Com'olutioii integral i)araiiioters 








x± 


1.21 ±0.06 


(2.82 ±0.20) 
+ (0.81±0.23)i 


(2.32 ±0.16) 

+ (0.50±0.18)i 


(3.44 ± 0.47) X 
- (3.91 ± 1.08) 



Quark and ly-boson masses 
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Mw 


1.2 ± 0.1 GoV 


1.2 ± 0.2 GoV 


0.10 ± 0.03 GoV 


172.5 ±2.7 GoV 


80. 1 GoV 



Table 1: The numerical values of the parameters. 
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where h{u, x) is the hard-scattering function: 
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and where Li2 is the usual dilogarithm function Li2(x) 
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$x is the hght-cone wave function with transverse polarization and $^1 is a distribution am- 
plitude of the B meson involved in the leading-twist projection. We have also: 
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The first negative moment of $bi can be parametrized by a quantity such as 



Following [1], and defining 
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= -6 a_i+ < v'^ >_L . 
The numerical values of all the needed parameters can be found in Table [H 
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Using the equations mentioned in this section, Superlso is able to calculate the isospin asym- 
metry of i? — >■ K*'^. For the computation of the inclusive branching ratio, we used the relations 
of pis]. 



3 Content of the Superlso package 

Superlso is a C program respecting the C99 standard (through the use of the complex. h 
library), whose main purpose is to calculate the isospin symmetry breaking of i? — > K*'-y decays 
in the MSSM with minimal flavor violation. In this package, five main programs are provided, 
but the users are also invited to write their own main programs. main_example . c is intended 
as an example for writing new programs, and slha. c can scan files written following the SUSY 
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Les Houches Accord (SLHA) [6J, and calculates the corresponding isospin asymmetry. The 
main programs msugra.c, amsb.c and gmsb.c have to be linked to the ISASUGRA/ISAJET [7J 
and/or the SOFTSUSY |E] packages, in order to compute mass spectra and couplings within the 
minimal Supergravity scenario (mSUGRA), the minimal Anomaly Mediated SUSY Breaking 
scenario (AMSB) or the minimal Gauge Mediated SUSY Breaking scenario (GMSB). 
The calculation of isospin asymmetry requires the following steps: 

• Generation of a SLHA file with ISAJET or SOFTSUSY, or the SLHA file provided by the 
user, 

• Scan of the SLHA file, 

• Computation of the required Wilson coefficients, 

• Calculation of the isospin asymmetry. 

Superlso also computes the inclusive branching ratio of 6 — > 57, and the code is modular enough 
so that routines calculating other observables, such as the branching ratio of Bg —>■ /^^A*" or the 
anomalous magnetic moment of the muon, can be easily added. For this purpose, a procedure 
has to be added in src/ which can also benefit from the functions already implemented in the 
program (for example for the calculation of the Wilson coefficients and the scan of the SLHA 
files.) 

The code includes the definition of a structure in src/include .h: 

typedef struct parameters 
{ 

int model; /* mSUGRA = 1, GMSB = 2, AMSB = 3 */ 

int generator; /* ISAJET = 1, SOFTSUSY = 2 */ 

float Q; /* Qmax ; default = M.EWSB = sqrt (m_stopl*mstop2) */ 

float m0,ml_2,tan_beta,sign_mu,A0,mass_W; /* mSUGRA parameters */ 

float Lambda, Mmess,N5,cgrav,m3_2; /* AMSB, GMSB parameters */ 

float mass_Z,mass_b,mass_top_pole,mass_tau_pole; /* SM parameters */ 

float inv_alpha_em, alpha_s_MZ, alpha, Gfermi,GAUGE_Q; /* SM parameters */ 

float charg_Umix [3] [3] , charg_Vmix [3] [3] ,stop_mix[3] [3] ,sbot_mix[3] [3] 

, stau_mix [3] [3] ,neut_mix [5] [5] ,mass_neut [5] ; /* mass mixing matrices */ 

float Min , Ml.Min , M2_Min , M3_Min , At_Min , Ab.Min , Atau_Min , M2Hl_Min , M2H2_Min , 

mu_Min,M2A_Min,tb_Min,mA_Min; /* optional input parameters at scale Min */ 

float MeL.Min , MmuL_Min , MtauL.Min , MeR.Min , MmuR.Min , 

MtauR_Min; /* optional input parameters at scale Min */ 

f loat MqLl_Min , MqL2_Min , MqL3_Min , MuR.Min , McR.Min , MtR.Min , MdR.Min , MsR.Min , 
MbR_Min; /* optional input parameters at scale Min */ 
float N51,N52,N53,M2H1_Q,M2H2_Q; /* optional input 
parameters (N51...3: GMSB) */ 

float mass_d,mass_u,mass_s ,mass_c ,mass_t ,mass_e ,mass_nue,mass_mu,mass_num. 
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mass_tau,mass_nut ; /* SM masses */ 

float mass_gluon,mass_photon,mass_ZO; /* SM masses */ 

float mass_hO,mass_cHO ,mass_AO,mass_H,mass_dnl ,mass_upl ,mass_stl ,mass_chl , 
mass_bl ,inass_tl ; /* Higgs & superparticle masses */ 

float mass_el ,mass_nuel ,mass_mul ,mass_numl ,mass_taul ,mass_nutl ,mass_gluino , 
mass_chal ,mass_cha2 ; /* Higgs & superparticle masses */ 

float mass_dnr ,mass_upr ,mass_str ,mass_chr ,mass_b2,mass_t2,mass_er,mass_mur , 

mass_tau2; /* superparticle masses */ 

float mass_nuer ,mass_numr ,mass_nutr ,mass_graviton, 

mass_gravitino ; /* superparticle masses */ 

float gp,g2,g3,YU_Q,yut,YD_Q,yub,YE_Q,yutau; /* Yukawa couplings */ 
f loat HMIX.Q , mu_Q , tanb.GUT , Higgs.VEV , mA2_Q , MSOFT.Q , M1_Q , M2_Q , 
M3_Q; /* parameters at scale Q */ 

f loat MeL_Q , MmuL_Q , MtauL_Q , MeR_Q , MmuR.Q , MtauR.Q , MqLl.Q , MqL2_Q , 
MqL3_Q,MuR_Q,McR_Q,MtR_Q,MdR_Q,MsR_Q,MbR_Q; /* masses at scale Q */ 
float AU_Q , A_u , A_c , A_t , AD_Q , A_d , A_s , A_b , AE_Q , A_e , A_mu , 
A_tau; /* trilinear couplings */ 
} 

parameters ; 

This structure contains all the important parameters and is called by several functions in the 
program. We can now review the important routines of the code. 

• void Init_param (struct parameters* param) 

This function initializes the param structure, setting the parameters at 0, apart for the 
SM masses and the value of the strong coupling constant at the Z-boson mass, which 
receive the values given in [9]. Moreover, model is set to -1 as an indicator for a wrong 
reading of the SLHA file. After the reading step, a value remaining -1 indicates a problem 
with the SLHA file. 

• int Les_Houches_Reader (char name [] , struct parameters* param) 

This routine is able to read a SLHA file whose name is contained in name, and put 
all the read parameters in the structure param. This function has been written based 
on It can be reused independently from Superlso, provided struct parameters is 
defined. 

• int sof tsusy_sugra(f loat mO, float ml2, float tanb, float AO, 
float sgnmu, float mtop, float mbot, float alphas_mz, char name [] ) 

• int isajet_sugra(f loat mO, float ml2, float tanb, float AO, 
float sgnmu, float mtop, char name[]) 
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int softsusy_gmsb (float Lambda, float Mmess, float tanb, int N5, 
float cGrav, float sgnmu, float mtop, float mbot, float alphas_mz, 
char name [] ) 

int softsusy_amsb (float mO, float m32, float tanb, float sgnmu, 
float mtop, float mbot, float alphas_mz, char name [] ) 

The above four routines call SOFTSUSY or ISAJET to compute the mass spectrum cor- 
responding to the input parameters. The input parameters are described in the next 
section. These routines return a SLHA file whose name has to be specified in the string 
name. 

float alpha_s_running (float Q, float mtop, float mbot, 
float alphas.MZ, float MZ) 

This function computes the strong coupling constant at the energy scale Q, provided the 
top quark mass mtop, bottom quark mass mbot, Z-boson mass (Mz) MZ and as{Mz) 
alphas_MZ are given, so that the matching can be made between the scales correspond- 
ing to different flavor numbers. The main formula for calculating is based on [9j. 

float running_mass (float quark_mass, float Qinit, float Qfin, 
float mtop, float mbot, float alphas_MZ, float MZ) 

This function calculates the running quark mass at the scale Qfin, for a quark of mass 
quark_mass at the scale Qinit, knowing the top quark mass mtop, bottom quark mass 
mbot, Z-boson mass (Mz) MZ and as{Mz) alphas_MZ, based on [TO] . 

void C.calculator (float C0[], float Cl[], float mu, 
struct parameters* param) 

This procedure computes the LO contributions to the Wilson coefficients Ci...s C0[] as 
well as the NLO contributions CI [] at the energy scale mu, using the parameters of param, 
and based on [T]. 

float deltaOm(float C0[], float CO.spec [], float CI [] , 

float Cl_spec [], struct parameters* param, float mub, float muspec, 

float lambda_h) 

This function computes the isospin asymmetry as described in the precedent section, 
using both the LO and NLO parts of the Wilson coefficients at scale mub (CO [] and 
Cl[]), and at the spectator scale muspec (CO_spec[] and Cl_spec[]), with the addi- 
tional input A/i lambda_h, in GeV. It uses the complex. h library which is available in C 



8 



Particle 




X? 


Ir 




xt 


ti 


9 


6i 


n 


Qr 


Lower bound 


111 


46 


88 


43.7 


67.7 


92.6 


195 


89 


81.9 


250 



Table 2: Lower bounds on the particle masses in GeV [9j. 
compilers respecting the C99 standard. 



• float deltaO_calculator (char name[]) 

This function is somehow a "container" function whose argument is the name of the 
SLHA file, and which calls first Init_param and Les_Houches_Reader, then calculates 
the Wilson coefficients with C_calculator, and finally returns the results of deltaOm. 
main_example . c provides an example of computing the isospin asymmetry without call- 
ing this function, but by calling deltaOm directly. This is useful if the user wants to 
provide his/her own Wilson coefficients as input. 



• float BRbsgainma_calculator (char name[]) 

This function, similar to the precedent one, computes the inclusive branching ratio of 
b ^ sj (based on |5j). It has to call two additional functions, Cem_calculator and 
BR_bsgaimna. 



• int charged_LSP_calculator (char name[]) 

This routine, whose argument is the name of the SLHA file, checks whether the lightest 
supersymmetric particle (LSP) is charged, and if so returns 0, and if not returns 1 (getting 
-1 means that a problem occurred). Note that if gravitino is the LSP, then this routine 
verifies the charge of the NLSP. 



• int excluded_mass_calculator (char name [] ) 

This routine, whose argument is the name of the SLHA file, checks whether the parameter 
space point is excluded by the collider constraints on the particle masses summarized in 
tableEl and if so returns 1, otherwise returns (getting -1 means that a problem occurred). 
The constraints in table [2] can be easily updated by modifying src/masses . c. 
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4 Compilation and installation instructions 



The Superlso package can be downloaded at 

: / /www3 . tsl . uu . se/ ~nazila/ superisoT^ 
After unpacking the package, the main directory, 
superiso_vX.X 

is created. It contains the src/ directory, in which all the source files can be found. The main 
directory contains also a Makefile, a README, five sample main programs (main_example . c, 
msugra. c, amsb . c, gmsb . c and slha. c) and one example of input file in the SUSY Les Houches 
Accord format (example . lha). The compiler options should be defined in the Makefile, as 
well as the path to the ISAJET isasugra.x and SOFTSUSY softpoint.x executable files, 
when needed. This package has been written for a C compiler respecting the C99 standard. In 
particular, it has been tested successfully with the GNU C Compiler and the Intel C Compiler 
on Linux 32-bits or 64-bits machines, and with SOFTSUSY 2.0.14 and ISAJET 7.75. Additional 
information can be found in the README file. 
To compile the library, type 
make 

This creates the libisospin. a in src/. Then, to compile one of the five programs provided in 
the main directory, type: make name 

where name can be main_example, msugra, amsb, gmsb or slha. This creates an executable with 
the .X extension. Note that slha does not need ISAJET or SOFTSUSY programs, but msugra, 
amsb and gmsb need at least one of them. Beside, main_example is somehow a test-program 
in which ISAJET and SOFTSUSY have been disabled by default, but their use can be restored 
by uncommenting #def ine USE.ISAJET and #def ine USE.SOFTSUSY at the beginning of the 
program. 

main_example .X has to be run without any argument, and performs the calculation of test- 
points using the main routines of the library. It can be used as a basis for writting the user's 
own main program. 

slha.x calculates the isospin asjTumetry and the inclusive branching ratio using the parame- 
ters contained in the SLHA file whose name has to be passed as input parameter. 

amsb . X, gmsb . x and msugra . x also calculate the isospin asymmetry and the inclusive branch- 
ing ratio, starting first by calculating the mass spectrum and couplings thanks to ISAJET (for 
msugra. X only) and/or SOFTSUSY within respectively the AMSB, GMSB and mSUGRA param- 
eter spaces. For all these programs, arguments referring to the usual input parameters have to 
be passed to the program. If not, a message will describe which parameters have to be specified. 

A more detailled presentation of the input and output of these modules is given in the next 
section. 
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5 Sample input and output 



The inputs and outputs of the five provided main programs are described in the following. 
However, a 0- value for the calculation of isospin asymmetry or branching ratio corresponds to a 
problem in reading the SLHA file, which either has not been generated by ISAJET or SOFTSUSY 
for any reason, or is not conformed to the SLHA format. 

• main_example . X 

This program is an example demonstrating the capabilities of Super I so. No input argu- 
ment is needed as the calls to the procedure are hardcoded. If 
#define USE_ISAJET 
#define USE.SOFTSUSY 

is uncommented at the beginning of the program, it will call first the functions isajet_sugra, 
softsusy_sugra, sof tsusy_amsb and softsusy_gmsb, with an adequate set of input 
parameters and calculate the isospin asymmetry and branching ratio thanks to the rou- 
tines deltaO_calculator and BRbsgaimna.calculator. It will then calculate the isospin 
asymmetry and branching ratio corresponding to the parameters of the example . lha file. 
Finally, input parameters {i.e. the Wilson coefficients, the masses of strange, charm, bot- 
tom and top quarks, the Z-boson mass Mz and ^^(M^)) corresponding to the Standard 
Model case are given directly in the program, and the isospin asymmetry is calculated 
using the deltaOm function. This last calculation is provided as an example of calcula- 
tions of isospin asymmetry without a SLHA file or without ISAJET or SOFTSUSY. If the 
#def ine's are uncommented, the output is for example 



deltaO.isaj et_sugra=0 . 101417 
BR.isaj et_sugra=0 . 000257 

deltaO.sof tsusy_sugra=0 . 100604 
BR_sof tsusy_sugra=0 . 000260 

deltaO_sof tsusy_amsb=0 . 096156 
BR_sof tsusy_ainsb=0 . 000286 

deltaO.sof tsusy_gmsb=0 . 079418 
BR_sof tsusy_gmsb=0 . 000393 

deltaO_example_slha=0 . 100568 
BR_exajnple_slha=0 . 000261 

deltaO_direct_SM=0 . 086136 
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• slha.x 

This program calculates the isospin asymmetry and inclusive branching ratio while read- 
ing the needed parameters in SLHA files. It also checks whether the SUSY parameter 
space point is excluded by the experimental mass limits. For example, the command 

./slha.x example. lha 

returns 

deltaO=0 . 100568 
BR=0. 000261 
excluded_mass=0 

• msugra.x 

This program generates the isospin asymmetry and inclusive branching ratio correspond- 
ing to the mSUGRA parameters generated by ISA JET and/or SOFTSUSY. It also checks 
whether the SUSY parameter space point is excluded by the experimental mass limits, or 
whether the LSP is charged. If only one of these generators is available the corresponding 
#def ine has to be commented. The necessary arguments to this program are: 

— mo: universal scalar mass at GUT scale, 

— mi/2: universal gaugino mass at GUT scale, 

— Aq: trilinear soft breaking parameter at GUT scale, 

— tan/3: ratio of the two Higgs vacuum expectation values. 

Optional arguments can also be given: 

— sign{fj,): sign of Higgsino mass term, positive by default 

— mf'*^: top quark pole mass, by default 172.5 GeV, 

— ffTbiWi^): scale independent b-quark mass, by default 4.2 GeV (option only available 
for SOFTSUSY), 

— as{Mz): strong coupling constant at the Z-boson mass, by default 0.1172 (option 
only available for SOFTSUSY). 
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With SOFTSUSY 2.0.14 and ISAJET 7.75, running the program with: 
./msugra.x 500 500 50 



returns 



deltaO.sof tsusy=0 . 100604 
BR.sof tsusy=0 . 000260 
charged_LSP_sof tsusy=0 
excluded_masses_sof tsusy=0 

deltaO_isajet=0 . 101457 
BR_isajet=0. 000257 
charged_LSP_isaj et=0 
excluded_inasses_isaj et=0 

• amsb.x 

This program generates the isospin asymmetry and inclusive branching ratio correspond- 
ing to the AMSB parameters generated by SOFTSUSY. It also checks whether the SUSY 
parameter space point is excluded by the experimental mass limits. The necessary argu- 
ments to this program are: 

— mo: universal scalar mass at GUT scale, 
~ "^3/2: gravitino mass at GUT scale, 

— tan (5: ratio of the two Higgs vacuum expectation values. 

Optional arguments are the same as for the mSUGRA case. With SOFTSUSY 2.0.14, run- 
ning the program with: 

./amsb.x 400 10000 30 -1 



returns 
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deltaO=0. 110040 
BR=0. 000224 
excluded_mass=l 

• gmsb.x 

This program generates the isospin asymmetry and inclusive branching ratio correspond- 
ing to the GMSB parameters generated by SOFTSUSY. The necessary arguments to this 
program are: 

— A: scale of the SUSY breaking in GeV (usually 10,000 - 100,000 GeV), 

— Mmess- messenger mass scale (> A), 

— N^: equivalent number of 5 + 5 messenger fields, 

— tan/3: ratio of the two Higgs vacuum expectation values. 

Optional arguments are the same as for the mSUGRA case, with an additional one: 

~ Ccrav (> 1): ratio of the gravitino mass, to its value for a breaking scale of A, 1 by 
default. 

With SOFTSUSY 2.0.14, running the program with: 
./gmsb.x 60000 80000 5 50 

returns 



deltaO=0. 079901 
BR=0 . 000390 
excluded_mass=0 



6 Results 

The Super I so code has been extensively tested, and compared to the calculation of the isospin 
asymmetry available in the literature for the Standard Model [2], and with the Wilson coeffi- 
cients and inclusive branching ratios calculated for example by micrOMEGAS [TT], and they 
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globally show a good agreement. 

Figure [1] shows the dependence of the isospin asymmetry in function of the parameters of 
the mSUGRA and AMSB parameter spaces, and reveals that isospin symmetry breaking can 
strongly constrain the supersymmetric parameters space, especially for high tan (3 value. 
A more complete study of the constraints by isospin asymmetry in different supersymmetric 
parameter spaces is presented in [T2] . 



7 Conclusion 

The Superlso package allows to calculate the isospin asymmetry of B K*'y in the MSSM 
with minimal flavor violation, for any file in the SLHA format. It can also call ISAJET or 
SOFTSUSY to generate a SLHA parameter file in the mSUGRA, AMSB or GMSB supersymmet- 
ric parameter space. The package also provides a routine to calculate the inclusive branching 
ratio of B ^ Xs'j to enable a comparison between the constraints on the supersymmetric pa- 
rameters from the branching ratio and from the isospin asymmetry. 
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